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1.0  Introduction 


Synthetic  Aperture  Radar  (SAR)  is  an  airborne  radar  capable  of  generating  high  resolution  images  of  surface 
target  areas  and  terrain.  The  term  Synthetic  Aperture  is  used  because  the  radar  utilizes  the  motion  of  the  aircraft, 
or  other  platform,  to  synthesize  the  effect  of  a  large  aperture  antenna  from  a  physically  small  aperture  antenna. 
As  a  result,  high  resolution  images  can  be  obtained  from  a  SAR  using  an  antenna  which  would  otherwise  produce 
poor  resolution. 

The  high  resolution  image  data  produced  by  a  SAR  is  often  employed  with  a  monopulse  radar  to  locate  fixed 
ground  based  target.  However,  as  with  any  measuring  device,  there  exists  some  "uncertainty"  or  "error"  in  target 
location.  By  combining  the  SAR  measurements  and  monopulse  radar  measurements  with  measurements  from 
other  sensors,  we  hope  to  reduce  the  level  of  uncertainty  in  target  location.  The  goal  here  is  to  bring  together 
information  from  SAR  and  navigation  sensors  such  as  global  positioning  system  (GPS)  receiver  and/or  an  inertial 
navigation  system  (INS).  This  report  focuses  on  the  integration  of  measurements  from  a  SAR  and  an  INS  in  a 
Kalman  filter. 

The  contents  of  this  report  may  be  summarized  as  follows: 

1)  describe  the  operation  of  a  SAR 

2)  describe  the  operation  of  monopulse  radar, 

3)  describe  how  a  target’s  location  is  computed  from  SAR  and  monopulse  radar  measurements, 

4)  present  error  models  for  both  SAR  and  monopulse  radar  measurements, 

5)  derive  the  Kalman  filter  measurement  error  models 

Many  of  the  points  developed  here  have  been  presented  in  other  documents  [1,  2,  3,  4,  5,  6,  7,  8].  The  intent  of 
this  document  is  to  bring  together  some  of  these  ideas  to  implement  SAR  Measurements  in  a  centralized  Kalman 
filter.  The  goal  is  to  enhance  the  overall  system  accuracy. 

A  thorough  understanding  of  how  SAR  measurements  are  obtained  is  needed  to  employ  SAR  measurements  in  a 
Kalman  filter.  Therefore,  the  first  few  sections  of  this  report  provide  the  basic  background  and  develop  models  for 
SAR  measurements.  In  Section  2,  a  brief  overview  of  SAR  processing  techniques  is  presented.  This  overview 
includes  a  description  of  unfocused  line-by-line  processing,  focused  line-by-line,  and  Doppler  processing. 
Highlighted  in  this  discussion  are  the  measurements  that  are  inherently  available  from  a  SAR  for  computing  a 
target's  location.  Also,  explained  are  the  differences  between  stripmap  mode  and  spotlight  mode  SAR.  Section  3 
includes  a  brief  discussion  on  monopulse  radar.  Section  4  provides  an  introduction  to  the  geometry  of  SAR  and 
monopulse  radar  measurements.  Presented  here  are  the  coordinate  frames  and  coordinate  frame  transformation 
used  throughout  this  report.  Also,  presented  in  this  section  are  some  basic  mathematical  relationships  needed  in 
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later  sections  on  measurement  models.  Presented  in  Section  5  is  a  method  for  computing  a  target's  location  from 
SAR  measurements  and  monopulse  radar  measurements.  Section  6  presents  a  truth  model  for  both  the  SAR  range 
and  range  rate  measurements  and  the  monopulse  radar  azimuth  and  elevation  measurements.  This  section 
includes  a  detailed  stochastic  model  for  all  errors  present  in  radar  measurements. 

To  generate  a  measurement  for  a  Kalman  filter,  the  INS  outputs  must  be  related  geometrically  to  the  SAR  and 
monopulse  radar  measurements.  Section  7  describes  how  to  compute,  from  INS  indicated  outputs,  quantities  that 
are  similar  to  each  radar  measurement.  These  computed  quantities  are  referenced  in  this  report  as  the  INS 
indicated  measurements.  In  essence,  this  section  describes  the  geometric  relationships  between  INS  outputs  and 
radar  measurements.  Section  8,  derives  the  linearized  INS  indicated  measurements.  This  includes  a  perturbation 
analysis  which  describes  the  sensitivity  of  the  INS  output  errors  on  the  computed  INS  indicated  measurements. 
The  linearized  measurement  equations  derived  in  Section  8  are  used  in  Section  9  to  derive  the  linearized 
measurement  equations  needed  in  a  Kalman  filter. 

Finally,  the  last  few  sections  provide  a  summary  of  the  radar  and  INS  error  models.  Section  10  summarizes  the 
dynamic  equations  for  the  INS  and  radar  error  states.  Section  1 1  summarizes  the  measurement  equations  needed 
to  implement  a  Kalman  filter.  Section  12  summarizes  a  typical  error  budget  for  a  SAR/INS  integration.  Section 
13  provides  some  final  remarks  on  many  of  the  points  derived  in  this  report. 


2 


2.0  Synthetic  Aperture  Radar  Background 

The  primary  function  of  a  SAR  is  to  generate  high  resolution  radar  images  of  ground  terrain.  These  SAR  images 
are  often  clear  enough  for  an  operator  to  recognize  man-made  targets  on  the  ground.  As  a  result,  SAR  images  are 
employed  to  achieve  the  first  task  in  targeting,  namely  target  designation.  The  operator  designates  the  target 
manually  placing  a  cursor  over  the  desired  target  in  the  image.  Once  the  target  is  designated,  SAR  systems  often 
employ  conventional  radar  techniques  to  actually  locate  the  target.  For  instance,  the  range  (r)  and  range  rate  (r) 
to  a  given  target  can  be  computed  by  the  following  equations: 
cT 


r  = ' 


(1) 


r  = 


cIa 

2F 


(2) 


where  c  is  the  speed  of  light,  T  is  the  round-trip  transit  time  of  the  transmitted  radar  pulse,  Fd  is  the  Doppler 
frequency  shift,  and  F  is  the  frequency  of  the  transmitted  signal. 


The  basic  SAR  geometry  is  shown  in  Figure  1.  The  SAR  is  characterized  by  a  wide  beam  antenna  which 
illuminates  a  large  area  on  the  ground.  For  what  is  assumed  an  instant  in  time,  the  SAR  transmits  a  radar  pulse 
and  samples  the  magnitude  and  phase  of  the  return  signal.  Since  radar  waves  propagate  at  nearly  a  constant 
speed  in  the  Earth's  atmosphere,  the  earliest  samples  correspond  to  the  points  on  the  ground  nearest  the  aircraft. 
Likewise,  the  return  from  more  distant  points  are  represented  by  later  samples.  The  samples  are  stored  in  vectors 
of  data  referred  to  as  range  bins.  Figure  1  illustrates  the  area  on  the  ground  for  which  magnitude  and  phase 
information  are  collected  by  the  SAR  for  a  given  transmitted  pulse. 
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As  the  aircraft  travels  along  a  given  ground  track,  the  process  of  illuminating  the  ground  with  a  radar  signal  and 
sampling  the  magnitude  and  phase  of  the  return  is  repeated  every  T  seconds  for  consecutively  overlapping  areas 
on  the  ground.  As  shown  in  Figure  2,  the  result  is  a  two  dimensional  array  of  data  which  contains  magnitude  and 
phase  information  of  the  radar  returns  in  both  range  and  azimuth  directions. 


It  is  possible  to  generate  a  radar  image  directly  from  the  two  dimensional  array  of  magnitude  and  phase  data. 
However,  due  to  the  use  of  a  very  wide  beamwidth  antenna,  any  given  magnitude  and  phase  sample  represents  the 
sum  return  of  all  illuminated  "points"  on  the  ground  of  equal  distance  from  the  aircraft.  Consequently,  the 
resulting  image  will  have  very  low  resolution.  SAR  data  processing  must  be  performed  to  generate  a  high 
resolution  SAR  image. 


Range  Bins 
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pulse  was  transmitted 
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Figure  2.  SAR  Data  Storage. 


One  method  of  SAR  processing  is  called  "unfocused  line-by-line  processing"  (often  call  "unfocused  SAR").  In 
unfocused  SAR,  the  radar  returns  are  thought  of  as  a  line  of  elemental  radiators  along  the  flight  path.  To 
synthesize  a  long  array  of  antennas,  a  vector  sum  of  the  elements  in  each  range  bins  is  generated.  This  operation 
corresponds  to  the  summing  performed  by  the  feed  structure  that  interconnects  the  radiating  elements  of  a  real 
array. 
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Figure  3  illustrates  the  signal  processing  required  for  unfocused  SAR.  Note  in  Figure  3,  that  the  magnitude  of  the 
vector  sum  is  generated  and  then  shifted  onto  the  display  memory.  This  process  is  repeated  for  every  new  radar 
pulse  (or  every  Ith  pulse  for  some  systems),  thus  generating  a  real  time  high  resolution  image. 

One  major  drawback  of  "unfocused  SAR"  is  that  it  does  not  perform  well  when  the  image  is  generated  at  close 
ranges  to  the  target  area.  In  unfocused  SAR,  the  length  of  the  synthesized  array  of  antennas  must  be  short  enough 
in  relation  to  the  range  of  the  swath  being  mapped  that  the  line  of  sight  from  any  one  point  at  the  swath’s  range  to 
the  individual  array  elements  are  essentially  parallel.  As  illustrated  in  Figure  4,  if  the  "synthetic  array  length"  is  a 
significant  fraction  of  the  swath's  range,  then  the  range  to  a  given  point  on  the  "synthetic  array  boresight" 
becomes  different  for  each  array  element.  These  differences  in  range  can  result  in  considerable  differences  in  the 
phase  of  the  returns  thus  limiting  the  performance  of  the  SAR.  While  the  azimuth  resolution  can  initially  be 
reduced  by  increasing  the  length  of  the  synthetic  array,  a  point  is  soon  reached  which  any  further  increase  in 
length  degrades  the  performance  of  the  SAR. 


R,  R,  R4  K 


Figure  3.  Unfocused  line-by-line  SAR  data  processing. 
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Direction  of  Flight 


Figure  4.  Distance  variations  for  unfocused  line-by-line  processed  SAR. 


The  limitation  on  array  length  of  the  unfocused  SAR  may  be  largely  removed  by  "focusing"  the  synthetic  array. 
Processing  of  this  type  is  called  "focused  line-by-line  processing"  (often  called  "focused  SAR”).  In  principle, 
focusing  is  achieved  by  simply  applying  the  appropriate  phase  correction  to  the  returns.  For  example,  the 
necessary  phase  correction  0k,  for  the  k*  radar  pulse  can  be  approximated  as: 

6k  =  y(2ARk)  (3) 

where  X  is  the  wavelength  of  the  transmitted  pulse  and  ARk  is  the  range  difference  of  the  k*  radar  pulse  to  a  point 
on  the  boresight  when  compared  to  the  synthesized  array  boresight  range.  The  (2ARk  term  is  used  in  Equation  (3) 
to  account  for  the  round-trip  travel  of  the  radar  pulse.  Notice  in  Figure  3,  that  the  following  relationships  exist: 


(R  +  ARk)2  =  R2  +  dk 

(4) 

R2  +  2RARk  +  (ARk)2  =  R2  +  dk 

(5) 

ARk  2 

(6) 

2RARk(l+-2p-)  =  dk 

dk 

ARk  =  for  ARk  «  2R 


(7) 


where  R  is  the  synthetic  boresight  range  and  dk  is  the  distance  of  the  k*  pulse  from  the  synthesized  boresight. 
Substituting  Equation  (7)  into  Equation  (3)  yields  the  following  approximation  for  the  phase  correction: 


9k~  %  R  ' 


(8) 
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Figure  5  illustrates  the  required  processing  for  focused  SAR.  As  the  return  from  a  transmitted  pulse  comes  in, 
they  are  stored  in  the  top  row  of  the  range  bins.  When  the  return  for  the  last  range  bin  has  been  received,  the 
contents  of  every  range  bin  is  shifted  down  one  row  and  the  bottom  range  bin  is  discarded.  Between  these  shifts, 
each  range  bin  is  read  serially,  appropriately  "focused"  (phase  shifted)  around  the  central  array  element,  and 
summed.  This  focusing  and  summing  process  is  often  called  "azimuth  compression".  As  was  performed  for 
unfocused  SAR,  a  magnitude  is  computed  for  the  azimuth  compression  result.  This  magnitude  is  then  shifted  into 
the  display  memory. 

One  major  drawback  of  unfocused  line-by-line  SAR  processing  is  its  computational  complexity.  Notice  that  every 
element  of  each  range  bin  must  be  phase  shifted  and  summed  to  generate  a  single  line  on  the  display.  Every  time 
the  radar  transmits  a  new  pulse,  it  must  perform  a  phase  correction  all  over  again.  If  the  synthesized  array  is  very 
long  (many  samples  in  a  range  bin),  then  computational  requirements  become  immense. 

The  computational  requirements  can  be  reduced  significantly  by  a  SAR  processing  technique  called  "Doppler 
Processing".  The  computational  burden  is  reduced  by  processing  data  in  parallel  for  many  lines  of  the  image 
rather  than  one  line  at  a  time.  In  Doppler  processing,  the  returns  from  different  azimuth  angles  are  isolated  with 
Doppler  filters. 


Figure  5.  Focused  line-by-line  SAR  processing. 
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To  explain  the  details  of  this  approach,  first  consider  how  Doppler  frequency  is  related  to  the  azimuth  direction. 
As  an  aircraft  approaches  a  point,  the  range  is  decreasing  and  the  range  rate  is  a  negative  value.  When  the 
aircraft  is  broadside  of  the  target,  the  range  rate  becomes  zero  and  then  increases  to  a  positive  value  as  the  aircraft 
begins  to  egress  from  the  point.  This  fact  is  illustrated  in  Figure  6. 


Figure  6.  Doppler  velocity  of  an  aircraft  as  it  passes  a  point  on  the  ground. 

The  Doppler  histories  of  several  evenly  spaced  points  at  the  same  offset  range  is  shown  in  Figure  7.  As  you  can 
see,  the  Doppler  histories  are  identical,  except  for  being  staggered  slightly  in  time.  Consequently,  at  any  point  in 
time,  the  return  from  every  point  has  a  slightly  different  frequency.  The  difference  in  frequency  corresponds  to  an 
azimuth  separation  of  the  points.  Therefore,  by  virtue  of  the  differences  in  the  Doppler  shift,  we  can  isolate  the 
return  received  from  each  point. 


-Ilvjf—  - - - - - - * - • 

Figure  7.  Doppler  history  of  evenly  spaced  points  on  the  ground. 
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Figure  8  illustrates,  in  block  diagram  form,  how  Doppler  processing  is  performed.  In  Doppler  processing,  a  phase 
correction  is  applied  at  the  onset  of  data  collection.  This  process,  called  focusing,  converts  the  return  from  each 
point  on  the  ground  to  a  constant  Doppler  frequency.  That  frequency  corresponds  to  the  Doppler  frequency,  as 
seen  from  the  center  of  the  segment  of  flight  over  which  radar  return  data  was  collected.  When  a  complete  set  of 
data  is  collected  in  all  range  bins,  the  data  is  applied  to  a  bank  of  Doppler  filters  which  separates  the  azimuth 
components  of  the  data.  The  magnitude  of  the  result  is  then  shifted  into  the  display  memory. 


At  the  end  of 
the  integration 
period,  outputs 
of  all  filters  arc 
transferred  to  the 
display  memory. 


Figure  8.  SAR  Doppler  processing. 


Notice  in  Doppler  processing  that  phase  corrections  are  applied  at  two  points  in  the  processing  algorithm,  namely, 
when  the  return  is  focused  and  when  the  Doppler  filtering  is  performed.  For  focusing,  only  one  phase  correction 
is  required  for  each  range  bin  element.  The  number  of  phase  corrections  needed  for  the  back  of  Doppler  filters  is 
on  the  order  of  Nlog2N,  where  N  is  the  number  of  integrated  radar  returns.  Thus,  the  total  number  of  phase 
corrections  required  for  Doppler  processing  is  on  the  order  of  N(1  +  log2N).  It  is  easy  to  recognize  that  this  value 

is  somewhat  less  than  the  N2  phase  corrections  needed  for  focused  line-by-line  SAR  processing.  Herein  lies  the 
major  advantage  of  Doppler  processing. 

The  SAR  processing  techniques  above  have  been  presented  with  the  assumption  of  stripmap  mode  SAR. 
Stripmap  SAR  is  simply  a  side-looking  radar  that  maps  a  strip  of  area  on  the  ground  parallel  to  the  flight  path. 
However,  by  gradually  changing  the  look  angle  of  the  real  antenna  as  the  aircraft  travels,  the  radar  can  repeatedly 
map  a  given  region  of  interest.  This  mode  of  operation,  called  spotlight  mode  SAR,  enables  the  operator  to 
maintain  surveillance  over  an  area  for  a  longer  period  of  time.  Also,  since  the  beam  dwells  continually  over  the 
mapped  area,  the  length  of  the  synthetic  array  is  not  limited  by  the  real  antenna  beamwidth.  Consequently, 
spotlight  SAR  is  capable  of  generating  superior  quality  images.  Figure  9  illustrates  the  differences  between 
stripmap  mode  SAR  and  spotlight  mode  SAR. 
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Figure  9.  Stripmap  mode  SAR  versus  spotlight  mode  SAR. 
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3.0  Monopulse  Radar  Background 

A  monopulse  radar  measures  the  azimuth  and  elevation  angles  of  the  line-of-sight  vector  from  the  aircraft  to  the 
target.  The  basic  vector  geometry  defining  the  azimuth  angle  (A)  and  elevation  (£)  angle  are  shown  in  Figure  10. 
Here,  the  azimuth  angle  is  defined  to  be  the  angle  between  a  vector  pointing  along  the  antenna  frame  y-axis  and 
the  line-of-sight  vector.  The  elevation  angle  is  defined  to  be  the  angle  between  a  vector  pointing  along  the 
antenna  frame  z-axis  and  the  line-of-sight  vector.  The  antenna  frame  is  defined  in  detail  in  the  next  section  of 
this  report. 


The  term  monopulse  refers  to  the  ability  of  the  radar  to  obtain  complete  angle  information  from  a  single  radar 
pulse.  In  a  monopulse  radar,  a  pair  of  antennas  is  used  to  form  two  similar  beams  which  point  in  slightly  different 
directions.  These  two  antenna  ports  are  connected  to  a  hybrid,  a  device  for  combining  the  receiver  RF  signals,  to 
obtain  a  sum  and  difference  output.  The  antenna  pattern  for  the  difference  port  will  exhibit  a  null  directly 
between  the  beams.  Likewise,  the  pattern  for  the  sum  port  will  exhibit  peak  directly  between  the  beams.  This 
makes  possible  a  pencil-beam  tracking  and/or  angle  measurement  radar  that  exhibits  very  high  precision. 

Many  SAR  systems  include  a  monopulse  radar  mode.  In  operation,  the  monopulse  processor  measures  the 
azimuth  and  elevation  angles  for  every  cell  in  the  SAR  image. 


11 


4.0  Measurement  Geometry 

In  this  section,  the  coordinate  frames  are  defined  and  the  coordinate  transformation  matrices  are  developed.  Also, 
this  section  summarizes  some  basic  geometric  relationships  for  SAR  range  and  range  rate  measurements  and 
monopulse  radar  azimuth  and  elevation  measurements. 

4.1  Coordinate  Frames 

The  coordinate  systems  of  interest  for  this  problem  are  the  Earth-centered  Earth-fixed  (ECEF)  frame,  the 
navigation  frame,  the  heading  frame,  the  path  frame,  the  body  frame,  and  the  antenna  frame.  These  right  handed 
orthogonal  coordinate  frames  are  defined  below.  Refer  to  Figure  1 1  for  a  graphical  interpretation  of  each  of  these 
frames. 

1)  Earth  centered  fixed  frame  (xe,  ye,  ze): 

As  the  name  suggests,  the  ECEF  frame  has  its  origins  at  the  Earth's  center  of  mass.  The  7?  axis  is  the  Earth's  axis 
of  rotation  and  points  toward  the  north  pole.  The  y?  -  ye  axis  lies  in  the  Earth's  equatorial  plane  such  that  x6  is 
fixed  at  Greenwich. 

2)  Navigation  frame  (xn,yn,  zn): 

The  navigation  frame  has  its  origins  at  the  aircraft’s  center  with  the  xn-yn  plane  tangent  to  the  reference  ellipsoids 
(i.e.  the  navigation  frame  is  a  local  level  frame).  In  this  frame,  xn  points  north,  >'n  points  east,  and  zn  points 
down  perpendicular  to  the  reference  ellipsoid.  The  location  of  the  navigation  frame  is  typically  specified  relative 
to  the  ECEF  frame  by  the  longitude  A,  the  latitude  <j>,  and  the  altitude  h. 

3)  Heading  frame  (xh,  yh,  zh): 

The  heading  frame  is  local  level  frame  whose  origin  is  located  at  the  aircraft's  center.  The  xh  axis  points  in  the 
direction  of  the  local  level  component  of  the  velocity  vector.  The  zh  axis  points  down  toward  the  reference 
ellipsoid  and  the  vh  is  tangent  to  the  reference  ellipsoid.  As  shown  in  Figure  11(b),  the  heading  frame  is  defined 
relative  to  the  navigation  frame  by  a  rotation,  denoted  J)z,  about  the  zh  axis. 

4)  Path  frame  (xP,  yP,  zP): 

The  path  frame  has  its  origins  at  the  aircraft's  center  such  that  the  xp  axis  points  in  the  direction  of  the  velocity 
vector.  The  yP  axis  is  tangent  to  the  reference  ellipsoid  and  points  in  the  same  direction  as  yh.  The  zP  point  in  an 
appropriate  direction  for  a  right-handed  orthogonal  coordinate  frame.  Note  in  Figure  11C,  the  path  is  defined 
relative  to  the  heading  frame  by  a  rotation,  denoted  rjy,  about  the  >,h  axis. 
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5)  Body  frame  (xb,yb,  zb): 

As  depicted  in  Figure  11D,  the  body  frame  defines  the  attitude  of  the  aircraft's  body  relative  to  the  navigation 
frame  via  roll,  and  yaw  angles  ax,  ay,  and  az.  Its  origin  is  at  the  aircraft's  center  such  that  xb  points  out  the  nose 

of  the  aircraft,  yb  points  out  the  right  side  of  the  aircraft,  and  zb  points  out  the  bottom  of  the  aircraft. 

6)  Antenna  frame  (xa,  ya,  za): 

The  antenna  frame  defines  the  orientation  of  the  radar's  antenna  relative  to  the  body  frame  via  antenna  roll  angles 
antenna  pitch  angle  Cy,  and  antenna  yaw  anglecr  The  antenna  frame  x3  axis  is  defined  such  that  it  points  in  a 

direction  perpendicular  to  the  antenna  phase  plane.  The  ya  axis  and  za  both  lie  in  the  phase  plane  of  the  antenna 
such  that  ya  points  along  the  semi-minor  axis  and  z2  points  along  the  semi-major  axis.  The  origin  of  the  antenna 
frame  is  assumed  to  be  at  the  aircraft's  center. 


4.2  Coordinate  Transformations 

Coordinate  transformations  between  frames  are  often  required  to  perform  vector  additions  and  multiplications. 
The  necessary  coordinate  transformations  are  summarized  below.  In  order  to  write  the  results  compactly,  we 
abbreviated  cos  by  C  and  sin  by  5. 


Navigation  frame  to  ECEF  frame 


-S<|>CA.  -SX  -C$CX  ' 
-S<(>SA,  CX  -C<j)SA. 
C<{>  0  -S<t>  . 


Heading  frame  to  navigation  frame 


Cftx 

Stix 

0 


-Sfix  o 

Ciix  o 
0  1  . 


Path  frame  to  heading  frame 


0 

1 

0 


Sfiy 

0 

CtV  . 


Body  frame  to  navigation  frame 


CazCay 

SazCccy 

-Say 


CazSaySax  -  SazCax  CazSayCax  +  SazSax 
SazSaySax  +  CazCax  SazSayCax  -  CazSax 
CaySax  CayCax 


(9) 


(10) 


(ID 


(12) 
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(13) 


Antenna  frame  to  body  frame 


C^C%- 

-S% 


^z^^y^^x  "  ^z^x 

S^S^S^X  +  C^Cqx 


+  S?zScx 

S^zS^yC?x-C?zS^ 

C<5yC(^ 


Figure  11.  Graphical  Depiction  of  the  Euler  angles  between,  A)  the  navigation  and  ECEF 
frames;  B)  the  heading  frame  and  the  navigation  frame;  C)  the  path  frame  and  the 
heading  frame;  D)  the  body  and  navigation  frames;  and  E)  the  antenna  and  body 
frames.  Note  in  (D)  and  (E),  the  frames  were  separated  for  clarity. 


4.3  Geometry  of  SAR  Range  and  Range  Rate  Measurements 

The  basic  vector  geometry  of  an  aircraft  and  a  ground  based  target  in  an  Earth-center  fixed  (ECEF)  coordinate  is 
shown  in  Figure  12. 
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g 

The  relationship  between  the  ECEF  referenced  line-of-sight  vector  pe,  the  target  location  vector  rt,  and  the 
6 

aircraft  position  vector,  ra,  is  given  by: 


P  e  e 
P6  =  rt  -  ra 


(14) 


The  range  from  the  aircraft  to  the  target  is  the  magnitude  of  p  which  is  found  by  computing  the  2-norm  of  the 
vector. 


r  =  II  pe  II  =  V[(Pe)T  Pel  (15) 


Figure  12.  Target  Location  Geometry. 


The  range  rate  is  found  by  computing  the  time  derivative  of  r. 
.  ^  dr  _  (Pe)T  ¥ 

r_dt~V[(pe)Tpe] 


(16) 


For  stationary  targets,  which  we  now  assume,  note  that  the  term,  in  Equation  (16)  becomes: 

dpf  _  drj  dr^_  e 
dt  -  dt  "  dt  ~  v® 


(17) 


where  va,  is  the  Earth-referenced  aircraft  velocity  expressed  in  the  Earth  frame.  Also,  note  that  the  unit  look 
vector  between  the  aircraft  and  the  ground  target  may  be  expressed  as: 

a  ne  ne 

6  _  P  _  P  (18) 

“p  ~  it  pe  11  -  Vt(pe)Tpe] 

Substituting  Equations  (17)  and  (18)  into  Equation  (16)  yields  the  following  expression  for  the  range  rate: 

f=-(uJ)Tvea  (19) 

It  is  easily  shown  that  Equations  (15)  and  (19)  can  be  expressed  directly  in  terms  of  the  navigation  frame  vectors, 
i.e.: 

r  =  II  pn  ll=V[(Pn)T  Pn]  (20> 
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(21) 


4.4  Geometry  of  Monopulse  Radar  Azimuth  and  Elevation  Angle  Measurements 

Recall  that  the  geometry  defining  the  azimuth  angle  (A)  and  elevation  (£)  angle  is  shown  in  Figure  10  above. 
Consider  the  unit  vectors  u^y  and  u^z  which  point  in  the  ya  and  za  directions,  respectively.  These  vectors  are 
expressed  in  the  antenna  coordinates  as 

i4  =  [0  10]T  (22) 

«i  =  [0  0  1]T  (23) 

and  may  be  converted  to  the  navigation  frame  by  employing  coordinate  transformations: 

%  =  CJC£>&  =  CJ,V  <24) 

^  =  =  «  (25) 

The  relationship  between  A  and  E  (quantities  measured  by  the  radar)  and  the  unit  vectors  of  interest  are  given  by: 
cos(A)=  (u£) T  u£y  (26) 

cos(E)=  (iy>T  ly  (27) 

Solving  for  A  and  E  yields 

A  =  cos"1  [(u£)T  Upy]  (28) 

E  =  cos-1[(^Sz]  (29) 

where  both  A  and  E  lie  in  the  interval  [0, 7t]  radians. 
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5.0  Targeting  with  SAR  Measurements 

In  general,  SAR  targeting  algorithms  need  three  pieces  of  information  to  compute  a  target's  location,  namely 
range,  range  rate,  and  height  above  the  target.  As  shown  in  Figure  13,  the  intersection  of  the  sphere  of  constant 
range,  cone  of  constant  range  rate  (centered  about  the  aircraft  velocity  vector),  and  the  plane  of  constant  height 
above  target  defines  the  location  of  the  target. 


Figure  13.  Range-Doppler  target  location  geometry. 


Range  and  range  rate  are  inherently  measured  and  stored  for  each  cell  in  a  Doppler  processed  SAR  image.  The 
height  above  target,  can  be  computed  from  the  range  measurement  and  the  precise  angle  measurements  of  a 
monopulse  radar.  Derived  next  are  the  mathematical  equations  needed  to  solve  the  target's  location  form  radar 
measurements.  For  the  development  which  follows  refer  to  Figure  14. 
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Here  we  consider  the  target’s  location  vector  in  path  frame  coordinates.  Notice  that  Doppler  cone  angle,  p,  is  the 
direction  cosine  of  the  line-of-sight  (LOS)  vector  relative  to  the  x-axis  of  the  path  frame.  Therefore,  the  following 
relationship  exists  between  the  Doppler  cone  angle  and  the  x  component  of  the  LOS  vector 


cos(p)=y 

(30) 

Solving  for  the  p*  yields 

Px  =  r  cos(p) 

(31) 

Note  that  if  p  <  n/2,  the  cos(/J)  can  be  computed  by  the  following  equation: 

IVtos^ 

COS(P)=  |  | 

ITa' 

(32) 

where  1  vLOgl  is  the  component  of  the  Earth  reference  velocity  va  along  the  radar  line-of-sight. 

Recognizing  the 

fact  that  1  vL0S  1  =  I  r  1  and  substituting  Equation  (32)  into  (31)  yields  the  following  expression  for  the  x- 

component  of  pP 

p-r!MI 

(33) 

The  y-component  of  pP  is  derived  using  the  trigonometric  identity,  sin2(P)  =  1  -  cos2(P)  ,  and  Equation  (30)  to 

obtain  the  following  expression: 

sin2(P)  =  l-^ 

(34) 

(35) 

<p')2  +  (Pz)2 
(Px)2  +  (Py)2  +  (Pz)2 

(36) 

_  (Py)2  +  (Pz)2 

sin2(p)  -  r2 

(37) 

Solving  for  (p^)2  : 

(Py)2  =  Vr2  Sin2(P)  -  (Pz)2 

(38) 

(Py)2  =  '\jr2  0  -  COS2(P))  -  (Pz)2 

(39) 
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In  Equations  (33)  and  (40),  we  now  have  expressions  for  the  path  frame  x  and  y  components  of  the  LOS  vector. 
The  known  quantities  in  these  equations  are  the  range,  r,  and  the  range  rate,  r,  (measured  by  the  SAR)  and  the 
aircraft  velocity  vector  va  (measured  by  the  inertial  navigation  system).  The  unknown  terms  are  the  components 

of  the  target  location  vector,  p£,  py,  and  pjj.  Hence,  we  have  two  equations  and  three  unknowns.  To  find  the  third 
equation  needed  to  solve  the  target's  location,  we  must  also  consider  the  heading  frame  coordinates. 


Note  in  Figure  14,  that  the  heading  frame  z-component  of  the  LOS  vector  is  equal  to  the  height  above  the  target 
measurement  such  that 

(41) 

Also,  recall  that  the  heading  frame  coordinates  are  related  to  the  path  frame  coordinates  by  a  simple  coordinate 
transformation,  i.e. 

(42) 

(43) 

(44) 


h  u 

Pz  =  ht 


Px  =  cos(ny)  Px  +  singly)  pi 

h  p 

Py  =  py 


pS  =  -sin(riy)  pi  +  cos(ny)  ppz 


From  Equations  (41)  and  (44)  it  is  easily  shown  that  the  path  frame  z-component  of  the  LOS  vector  is  given  by 
the  following  expressions: 


1  [  sin(r|  )  ppz  +  ht] 


Z  COS(T)y) 


Pz  =  tan(T|y)  r 


Iff 


lv. 


COS(Tly) 


(45) 

(46) 

(47) 


Using  the  results  derived  in  Equations  (33),  (40),  and  (47);  the  overall  vector  pP  is  expressed  as: 


Pp= 


Px 


L  Pz 


COS(T|y) 


(48) 


The  terms  cos(r|y)  and  sin(riy)  in  Equations  (48)  are  found  by  using  the  INS  indicated  velocities  in  the  following 


equations: 


cos(t)  )  = 


r  l  v„ 


-V, 


tan(ri  )  =  -  — 


(49) 


(50) 


where  vx  and  vz  are  the  heading  frame  x  and  y  components,  respectively,  of  the  aircraft  velocity  vector.  Next  we 
consider  the  computation  of  the  height  above  target,  ht,  from  monopulse  radar  angle  measurements. 


Given  the  monopulse  radar  measurements  of  azimuth  and  elevation  angle,  it  is  possible  to  compute  a  unit  LOS 
vector,  u^,  in  antenna  coordinates.  This  is  done  with  the  following  equations: 

u =  V 1  -  cos2(A)  -  cos2(E)  (51) 

u£y  =  cos(A)  (52) 

u^z  =  cos(E)  (53) 

where  =  [ujjx  u£y  u^z]T.  This  antenna  frame  unit  vector  may  be  expressed  in  navigation  frame  coordinates  by 
employing  the  appropriate  coordinate  transformations,  i.e. 

^=<^4  (54> 

Notice,  that  the  vector  u£  is  also  a  unit  vector.  Consequently,  the  component  of  Up  are  the  direction  cosines 
angles  (i.e.  the  cosine  of  the  angles  between  the  navigation  frame  coordinate  axes  and  the  line-of-sight  vector, 
denoted  cos(tp).  Using  this  fact  it  is  easily  shown  that  the  height  above  the  target  is  given  by: 

ht  =  r  cos(tp)  =  rujz  (55) 

where  uj)z  is  the  z-component  of  u"r 

As  shown  in  this  section,  it  is  possible  to  compute  a  target's  location  using  only  one  set  of  SAR  and  monopulse 
radar  measurements.  When  a  Kalman  filter  is  used  with  multiple-look  measurements,  it  may  be  possible  to  locate 
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a  target  with  a  smaller  set  of  measurements  (e.g.  range  measurements  only).  Also,  by  utilizing  the  complete 
measurement  history,  a  Kalman  filter,  provides  a  means  to  continually  improve  the  estimate  of  the  target's 
location. 
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6.0  Radar  Measurement  Models 


In  this  section,  we  present  the  measurement  equations  for  both  the  SAR  and  monopulse  radar  measurements. 
Here,  we  describe  the  different  sources  of  errors  present  in  each  measurement.  Also,  each  error  is  characterized 
by  a  Stochastic  model  so  that  they  may  be  incorporated  in  a  Kalman  filter  design. 

6.1  SAR  Measurements 

As  previously  mentioned,  a  SAR  is  capable  of  measuring  range  and  range  rate  to  a  given  target.  However,  these 
measurements  are  not  perfect  and  must  be  modeled  as  the  sum  of  true  value  and  the  measurement  errors  and 

noises.  The  measured  range,  denoted  r,  and  range  rate,  denoted  r,  are  modeled  by  the  following  equations: 


r  =  r  +  5rr.  +  r5C  -  v  -  v.  -  vr 

^  rDES  rQ  T 

r  =  ?  +  5?n  +  f  8C  +  ?8F  -  v*  -v* 

D  rDES  rQ 


(56) 

(57) 


where 

5rCL  is  the  range  clock  (phase)  error, 

5rD  is  a  Doppler  measurement  error, 

5C  is  the  radar  wave  propagation  speed  error 
5F  is  the  radar  frequency  error 

v  is  the  range  designation  error 

rDES 

v;  is  the  range  rate  designation  error 

rDES 

vr^  is  the  range  quantization  error 
v;  is  the  range  rate  quantization  error 
vr  is  the  range  timing  error 

rT 

vi  is  the  range  rate  timing  error 

rT 

V  is  range  error  which  captures  the  effects  of  multipath,  high  maneuvers,  and  edge  of  map  targeting 

rOT 

v*  is  range  rate  error  capturing  the  effects  of  multipath,  high  maneuvers,  and  edge  of  map  targeting. 

rOT 


Other  potential  range  rate  measurement  errors,  which  are  not  modeled  here,  include  the  errors  associated  with 
wave  movement  when  targeting  is  performed  over  bodies  of  water.  Presented  next  is  a  brief  description  of  each 
error  source. 
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Range  measurements  are  made  by  measuring  the  time  difference  between  the  transmitted  pulse  and  its  return. 
Consequently,  any  error  found  in  the  clock  rate  will  result  in  a  range  measurement  error.  In  general,  a  clock  runs 
either  slightly  fast  or  slightly  slow  with  very  slowly  changing  variations  in  the  rate.  Therefore,  the  error  8rCL  in 
the  range  measurement  resulting  from  a  clock  error  is  modeled  here  as  a  random  bias,  i.e. 

8?cl  =  0  (58) 

The  rms  value  of  this  error  is  on  the  order  of  0.01  feet. 

Recall  that  the  radar  measures  range  rate  by  measuring  the  Doppler  frequency  of  return  pulses.  However,  some 
radar  system  can  induce  Doppler  frequency  shifts  when  transmitting  radar  pulses  and  when  receiving  the  return. 
This  results  in  an  error  in  the  range  rate  measurement.  The  error  8rD  in  the  range  rate  measurement  resulting 
from  system  induced  Doppler  shift  error  is  modeled  here  as  a  random  bias,  i.e., 

8fD  =  0  (59) 

The  rms  value  of  this  error  is  on  the  order  of  0.001  feet/sec. 

Knowledge  of  the  radar  wave  propagation  speed  is  needed  to  map  the  measured  radar  wave  transit  time  to  a  range 
measurement  and  the  measured  Doppler  frequency  to  the  radar  range  rate  measurement.  Therefore,  as  shown  in 
Equations  (56)  and  (57),  the  propagation  speed  error  has  "scale  factor"  type  effect  on  both  the  range  and  the  range 
rate  measurements.  Likewise,  since  knowledge  of  the  radar  frequency  is  needed  to  map  Doppler  frequency  to  a 
range  rate  measurement,  error  in  the  radar  pulse  frequency  will  also  have  a  scale  factor  effect  on  the  range  rate 
measurement. 


The  propagation  speed  error,  8C,  and  frequency  error,  8F,  are  modeled  here  as  random  bias,  i.e. 

SC  =  0  (60) 

8F  =  0  (61) 

The  rms  value  of  the  propagation  speed  error  is  on  the  order  of  10  PPM.  The  rms  value  of  the  frequency  error  is 
on  the  order  of  20  PPM. 

The  target  designation  error  results  from  the  human  operator's  inability  to  designate  the  exact  pixel  in  the  SAR 
image  that  represents  the  desired  target.  The  target  designation  error  can  be  modeled  as  a  function  of  four 
parameters,  namely  clutter-to-noise  (CNR),  signal-to-noise  (SNR),  signal-to-clutter-plus  noise  ratio  (SCNR),  and 
the  resolution  cell  size.  The  CNR  establishes  that  the  image  quality  is  sufficient  for  individual  fixed  targets  to  be 
designated.  If  the  CNR  is  less  than  3  dB,  the  variance  of  the  target  designation  error  in  both  the  range  and 
azimuth  directions  is  assumed  large.  The  SCNR  is  compared  to  a  threshold  of  (lOdB)  to  establish  that  the 
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individual  targets  are  recognizable  by  a  human  operator.  The  CNR,  SNR,  and  SCNR  are  characterized  by  the 
following  expressions: 


CNR  =  pcAG2LaLrr1^  (62) 

SNR  =  PtAG2La Lr  Tx  (63) 

SCNR  =  CNR +  0.1  (64) 

where 

La  =  lO^216  x  10  =  atmospheric  loss  (65) 

Lf  =  10<5'4  x  10  6^r  =  rain  attention  (66) 

pc  =  0.1(Ar)(Aaz)  =  clutter  cross  -  section  (67) 


and  where  Tx  is  the  SAR  integration  time  AG2  is  the  two  way  loss  in  the  mainlobe  gain,  pv  is  the  target  cross- 
sectional  area,  7]  is  the  rain  rate  in  mm/hr,  r  is  the  range  to  the  target,  Ar  and  Aaz  are  the  azimuth  and  range 
resolutions,  respectively,  and  r0  is  the  range  at  which  the  SNR  is  unity  for  1  m2  target,  1  sec  array  time,  and 
undegraded  mainlobe  antenna  gain.  A  typical  value  of  r0  is  1.5  x  105  m .  For  a  stripmap  SAR,  range  and  azimuth 
resolutions  are  given  in  [9]  by  the  following  expressions: 


c  cT 
=  2B=  2 


(68) 


Aaz  =  f  (69) 

respectively,  where  c  is  the  speed  of  light,  B  is  the  bandwidth  of  the  transmitted  pulse,  T  is  the  unmodulated  pulse 
length,  and  D  is  the  antenna  diameter.  For  a  spotlight  mode  SAR,  the  range  resolution  is  the  same  as  above  and 
the  azimuth  resolution  is  expressed  as: 

=  i  (70) 

where  X  is  the  wavelength  of  the  transmitted  pulse,  and  <j)  is  the  angle  between  the  aircraft  and  the  target  created 
by  the  motion  of  the  aircraft  during  SAR  image  data  collection. 


In  earlier  discussion  SAR  processing,  it  was  shown  that  the  range  rate  measurement  is  directly  related  to  the 
azimuth  direction  of  the  SAR  image.  Consequently,  the  range  rate  resolution,  denoted  Ar,  is  a  function  of  the 
azimuth  resolution.  In  [9],  this  relationship  is  given  by  the  following  equation 

Ar  =  cotAaz  (71) 
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where  cot  is  the  angular  "rotation  rate"  of  the  aircraft  around  the  target.  The  angular  rotation  rate  is  given  by  the 
following  equation: 


lv"l 

wt=T 


i .HI 

'Va'j 


(72) 


When  the  (SCNR  >  10  dB),  the  designation  error  is  often  modeled  as  white  noise  with  a  range  component, 
denoted  v  ,  and  a  range  rate  components,  denoted  v*^.  The  rms  value  for  the  range  and  the  range  rate 


components  are  given  by: 

rms  value  for  v_  =  KAr 

rDES 

rms  value  for  v-  =  KAr 

J  rDES 


(73) 

(74) 


respectively,  where  typical  values  of  K  are  on  the  interval  [0.2,  2.5].  Generally,  K  is  inversely  proportional  to  the 
SCNR.  In  [2],  K  is  modeled  by  the  following  expression: 

1 


K  2y/SCNR 
Here,  we  will  model  K  as  a  constant  of  value  2.0. 


(75) 


The  range  quantization  error,  vf  ,  and  the  range  rate  quantization  errors,  v-^  are  a  result  of  range  and  range  rate 

resolution  errors.  In  [10],  quantization  error  is  modeled  as  a  uniformly  distributed  white  noise.  Therefore,  the 
range  and  range  rate  quantization  error  are  modeled  as  white  noises  with  an  rms  value  equal  to  the  following: 

Ar 


rms  value  for  vr^  =  ■ jj - 
rms  value  for  ^ 


(76) 


(77) 


The  range  timing  error,  v^,  and  the  range  rate  timing  error  v^,  result  from  the  time  delay  needed  to  perform  SAR 

processing.  These  timing  errors  are  modeled  as  white  noises  with  rms  values  given  by  the  following  expressions 
rms  value  for  vr  = Tdr  (78) 


rms  value  for  v-r  =  Tdr 


(79) 


respectively,  where  Td  is  the  delay  time.  The  time  delay  for  this  process  is  on  the  order  of  1  ms . 


The  white  noise  errors,  vr  ,  and  v;  ,  are  added  to  the  range  and  range  rate  measurements,  respectively,  to 

rOT  rOT 

capture  the  effects  of  other  miscellaneous  sources.  Such  errors  included  high  maneuver,  edge  of  map  targeting, 
signal  multipath,  and  small  squint  angles.  These  errors  are  not  always  present  in  the  SAR  measurements.  The 
lump  sum  effect  of  these  errors  are  modeled  here  as  white  noises  with  rms  values  given  by  the  following  equation: 
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rms  value  for  v  =  0,5 Ar 

(80) 

rms  value  for  v*  =  0.5 A  r 

(81) 

6.2  Monopulse  Radar  Measurements 

Recall  that  monopulse  radar  is  capable  of  measuring  the  azimuth  angle,  A,  and  the  elevation  angle  E.  Similar  to 
the  SAR  measurements,  the  monopulse  radar  measurements  are  also  corrupted  by  a  number  of  errors.  As  a  result, 
they  may  be  expressed  as  the  sum  of  the  true  value  and  the  measurement  errors. 


For  an  electronically  scanned  array  antenna  (ESA),  which  we  now  assume,  the  measured  azimuth  angle,  denoted 
A,  and  elevation  angle,  denoted  E,  are  modeled  by  the  following  equations: 


A  =  A  +  6At 


E  -  E  +  SEpemp  -  vEma  -  -  v£pu  -  v£cc  -  v£ms  -  v^  -  vEq  -  vEra 


where 

8ATEmp  and  $etemp  m  ^  temperature  induced  errors 

v A  and  Vr  are  the  mechanical  alignment  errors 
ama  c ma 

v  a  and  are  the  flexure  errors 

*FL  Z'FL 

v  a  and  vF  are  the  phase  uniformity  errors, 

A PU  CPU 

vA  and  Vr  are  the  Rx  channel  cross  coupling  terms 
*cc  a cc 

vA  and  vE  are  the  monopulse  slope  errors, 

MS  MS 

Va  and  Vr  are  the  thermal  noise  errors 
Am  c tn 

vA  and  vE  are  the  measurement  quantization  errors, 

vA  and  vF  are  the  random  deflection  errors. 
arad  crad 


Other  potential  errors  in  the  azimuth  and  elevation  angle  measurement,  which  are  not  modeled  here,  include 
atmospheric  bending  effects.  Presented  next  is  a  brief  description  of  each  of  the  above  errors. 


The  temperature  errors,  <54  jp^jp  and  SEj^p,  result  from  a  temperature  gradient  across  the  face  of  the  antenna 
array.  The  antenna  temperature  induced  errors  are  modeled  as  biases,  i.e. 


&atemp  ~  0 


(84) 


®  TEMP  “  0 


(85) 
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The  rms  value  for  the  temperature  induced  errors  are  on  the  order  of  0.4  mrads. 


To  increase  the  potential  scan  area,  many  SAR  antennas  are  designed  so  that  they  are  steered  both  electronically 
and  mechanically.  In  this  case,  there  will  be  pointing  errors  resulting  from  inaccuracies  in  both  the  electronic  and 
mechanical  steering  mechanisms.  The  errors,  vA  and  vE  ,  are  a  result  of  mechanically  steering  errors.  Other 

errors  in  this  section  capture  the  effects  of  electronics  steering  inaccuracies.  The  mechanical  alignment  errors  are 
modeled  here  as  white  noise  with  an  rms  value  on  the  order  of  0.3  mrads. 

The  flexure  errors  ,  vA  and  Vr  ,  model  the  effects  of  structure  vibrations  and  bending  between  the  antenna  and 

afl  cfl 

the  body  of  the  aircraft.  The  flexure  errors  are  modeled  as  white  noises  with  an  rms  value  on  the  order  of  0. 1 
mrads. 

The  positioning  errors,  vA  and  vEprf  are  due  t0  sources  such  as  phase  quantization  and  phase  shifter  errors. 
They  are  modeled  as  white  noises  with  an  rms  on  the  order  of  0.35  mrads. 

Cross  coupling  between  the  sum  and  difference  channel  of  the  monopulse  radar  can  result  in  a  pointing  error  of 
the  monopulse  beam.  The  errors  vA  and  vE  9  resulting  from  channel  cross  coupling  are  modeled  as  white 

noises  with  an  rms  value  on  the  order  of  0.15  mrads. 

Non-uniform  gains  across  the  face  of  the  antenna  array  can  also  result  in  radar  beam  pointing  errors.  These  errors 
are  often  referred  to  as  "monopulse  slope"  errors.  The  monopulse  slope  errors,  denoted  vA  and  vE  ,  are 

MS  MS 

modeled  as  white  noises  with  an  rms  value  on  the  order  of  0.1  mrads. 

In  [11],  thermal  noise  is  modeled  as  white  noise  process.  Therefore,  the  thermal  noise  errors,  denoted  vA  and 
vp  ,  are  modeled  here  as  white  noises  with  an  rms  value  of  0.2  mrads. 

CTN 

Quantization  errors  will  enter  into  the  azimuth  and  elevation  angle  measurements  as  a  result  of  the  finite  word 
length  used  in  the  digital  computers  which  process  the  monopulse  radar  measurements.  As  mentioned  earlier, 
quantization  errors  are  often  modeled  as  white  noise  processes.  Consequently,  the  azimuth  and  elevation  angle 
quantization  errors,  denoted  vA  and  vE^,  are  modeled  as  white  noises  with  an  rms  value  of  the  order  of  0.2  mrads. 

The  true  rms  value  for  this  error  is  highly  dependent  on  the  design. 
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The  antenna  random  will  often  bend  the  boresight  of  the  radar.  This  bending  is  usually  compensated  by  the  radar 

processor;  however,  often  there  exists  some  uncompensated  bending.  The  uncompensated  random  bending 

errors,  denoted  v,  and  vr  ,  are  modeled  as  white  noises  with  an  rms  value  on  the  order  of  0.4  mrads. 
arad  c rad 
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7.0  INS  Indicated  Range  and  Range  Rate 

The  INS  indicated  measurements  of  range  r  and  range  rate  r  may  be  expressed  by  the  following  relationships: 

r  =  V(Pn)T  Pn  (86) 

f=-(^)TVa  (87) 

where  pn  and  u”  are  computed  form  the  INS  indicated  aircraft  position  r"  and  the  indicated  target  position  r". 

p  "  =  (88) 


V  (pn)T  Pn 

Notice  that  an  indicated  target  position,  denoted  r”  is  needed  to  compute  the  INS  indicated  measurement  of  range 

and  range  rate.  However,  in  general,  it  is  not  provided  by  the  INS.  The  indicated  target  position,  r",  may  be 
obtained  using  a  pre-briefed  target  position.  If  the  pre-briefed  target  information  is  not  available,  simply  compute 
the  indicated  target  position  from  the  first  set  of  radar  measurements. 

INS  indicated  Azimuth  and  Elevation  Angle  Measurements 

The  indicated  values  of  the  unit  vectors  u^y  and  u^z  are  obtained  directly  form  Equations  (22)  and  (23),  i.e. 

u^y  =  (0  1  0]T  (90) 

u£2  =  [  0  0  1  ]T  (91) 

Likewise,  they  may  be  converted  to  the  navigation  frame  by  employing  coordinate  transformations  to  obtain: 

=  =  (92) 

<«) 

The  basic  relationship  between  INS  indicated  measurements  A  and  E  and  the  unit  vector  of  interest  is 

cos(A)=  (u£)T  u£y  (94) 

cos(E)=  (Up)T  iy  (95) 


Solving  for  A  and  E  yields: 
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A  =  cos'1  [(uJ)TuJy] 

(96) 

E  =  cos’1  [(u£)Tu£z] 

(97) 
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8.0  Linearized  INS  Indicated  Measurements 

Derived  in  this  section  are  the  linearized  equations  for  the  INS  indicated  measurements.  These  equations  will  be 
used  later  to  derive  the  Kalman  filter  measurement  equations.  Since  coordinate  transformations  are  required  to 
compute  INS  indicated  measurements,  this  section  begins  with  a  perturbation  analysis  of  the  INS  indicated 
transformation  matrices. 

8.1  Transformation  Matrix  Perturbations 

The  indicated  transformation  matrices  can  be  expressed  as  the  sum  of  the  true  transformation  matrix  and  the  INS 
indicated  transformation  matrix  error,  i.e. 


(fb  =  <fb+  6<fb 

(98) 

Q  =  c\+  8&z 

(99) 

Next,  we  find  an  expression  for  the  INS  indicated  body  to  navigation  frame  transformation  matrix  error  5c£.  A 

similar  development  applies  for  the  antenna  to  body  frame  transformation  error  <5c£.  Since  it  is  not  needed  for  the 
analysis  of  relative  SAR  targeting,  we  do  not  consider  the  perturbations  of  the  Earth  to  navigation  frame 
transformation,  heading  to  navigation  frame  transformation,  and  path  to  heading  frame  transformation. 

The  errors  in  the  INS  indicated  body  to  navigation  frame  transformation  matrix  are  a  result  of  errors  in  the  INS 
indicated  Euler  angles  (ax,  cya, )  where, 


ax  =  ax  +  5ccx 

(100) 

1 

ii 

+ 

(101) 

«2  =  az  +  fcz 

(102) 

Substituting  the  INS  indicated  Euler  angles  in  Equations  (100),  (101),  and  (102)  for  the  true  Euler  angles  in 
Equation  (12)  yields  the  following  expression  for  the  INS  indicated  transformation  matrix. 


C(az)C(ay)  C(az)S(ay)S(ax)  -  S(az)C(ax)  C(az)S(t^)C(ax)  +  S(az)S(ax) 
S(az)Cay  S(az)S(ay)S(ax)  +  C(az)C(ax)  S(az)S(^)C(ax)  -  C(az)S (ax) 
-Say  CayS(ax)  CayC  (ax) 


(103) 
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For  "small"  error  angles  (less  than  5  degrees),  the  following  approximations  exits  for  the  sin  and  cosine  of  an 
angle  plus  its  error. 

sin(9  +  80)  ~  sin(0)  +  80  cos(80)  (104) 

cos(0  +  80)  ~  cos(0)  -  80  sin(0)  (105) 

Using  equations  (104)  and  (105)  in  Equation  (103),  expanding  terms,  and  neglecting  and  2nd  and  3rd  order  error 
terms,  we  obtain  the  following  expression  for  the  INS  indicated  transformation  matrix  error. 


-5<*z 

+5azCazCay  -  Sc^Sc^Sc^ 


-do^Sa^OySc^  +  5  ctyCc^COySc^  +  So^Cc^SoyCc^  -  baJZo^Cc^  +  Sc^Sa^Sa^ 
+5azCazSaySax  +  5a^Sci^CaySc^  +  ba^Sa^a^Ca^  -  bo^Sc^Cc^  -  Sc^Co^So^ 
-  SOySOySa^  +  8  a^CayCc^ 


^azSazSayCccx  +  5ayCazCayC<\ '  KCazSccyS<h  +  8aZCazSoXCoX  +  HS«zCax 
+bazCazSayCax  +  bo^S^Cc^Cc^  -  Sa^So^Sc^Sc^  +  ba^o^So^  +  bo^Cc^Co^ 

-SOySGyCc^  -  Sa^COySc^ 


(106) 


Due  to  their  manipulative  properties  it  is  very  desirable  to  work  with  orthogonal  matrices.  However,  notice  that 
the  expression  for  SC£  in  Equation  (106)  is  nonorthogonal.  In  Britting  [12],  it  is  shown  that  the  "optimal" 

orthogonal  approximation  to  c£  in  the  sense  that  the  trace  of  [(C£)0  -  c£]T  [(CjJ)0  -  C^]  is  minimized  is  given  by: 


Substituting  Equations  (98)  into  Equation  (107)  and  expanding  the  square  root  term  in  series,  the  result  is 


(Cb)0  =  (I +  0.5[5C£  (C£)T  -  c£(5c£)T])  c£  (108) 

where  I  is  a  3x3  identity  matrix.  Since  it  involves  the  difference  between  a  matrix  and  its  transpose,  the  bracketed 
term  on  the  right-hand  side  of  Equation  (108)  is  skew  symmetric.  Consequently,  Equation  (108)  can  be  written 
as: 


(C$0=  {!-¥}< 


where 


»F= 


0 

V'z 

-Yy 


-V'z  V'y 

o  -v/x 

V'z  0  . 


(109) 


(110) 
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(111) 

(112) 
(113) 


The  non-zero  elements  of  'V  are  given  by: 

V'x  =  0-5[c315c21  -  c2j8c31  +  c328c22  -  c225 c32  +  c335c23  -  c238c33] 

V'y  =  0.5[cn8c31  -  c318cn  +  c138c32  -  c328c13  +  c138c33  -  c338c13] 

Yz  =  0.5[c218cjj  -  Cjj8c21  +  c225c12  -  c128c22  +  c235cI3  -  c138c23] 
where  c-  and  5cy  are  the  elements  of  cjj  and  Sc£,  respectively,  such  that  i  is  the  row  number  and/  is  the  column. 

Similar  to  the  above  development,  the  indicated  antenna  to  navigation  transformation  matrix  may  be  computer  by 
the  following  equation: 

(CJ)0={I-B}CJ  (114) 

where 

'  0  -Pz  Py  ' 

B=  Pz  0  -/?x  (115) 

.  -Py  Pz  0  . 

For  a  Kalman  filter  design,  it  is  sufficient  to  work  directly  with  the  error  angles,  yrx,  yfy,  yz,  Px,Py,  and  /3Z. 

However,  to  keep  in  mind  that  INS  systems  often  provides  only  the  Euler  angles.  Therefore,  it  may  be  necessary 
to  compute  the  Kalman  filter's  estimate  of  the  Euler  angels.  Therefore,  it  may  be  necessary  to  compute  the 
Kalman  filter's  estimate  of  the  Euler  angel  errors  from  the  Kalman  filter's  estimate  of  Vx,Yy,  and  \jfz,  to  perform  a 
correction  to  the  INS  indicated  Euler  angels.  For  example,  this  may  be  done  for  the  body  to  navigation  Euler 
angles  by  equating  Equation  (106)  to  'Pc£  and  solving  for  Sar  8ay,  and  8az,  and  in  terms  of  YrYy,  and  yz. 


8.2  Linearized  INS  Indicated  Range  and  Range  Rate  Measurements 

The  INS  indicated  position,  r",  and  the  indicated  target  position  r"  can  be  express  as  a  sum  of  the  true  values  and 
the  error  values,  i.e. 


—  n  n  c  n 

r,  =rt  +  5rt 

(116) 

ra  =  ra  +  5ra 

(117) 

Substituting,  Equations  (116)  and  (117)  into  Equation  (88)  yields 

^  n  c*n  n 

p  n  =  pn  +  gpn  =  ^  +  grt  .  ^  +  gra 

(118) 

where  5pn  =  8r"  -  Sr” 
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Express  the  INS  indicated  range  r  measurement  in  terms  of  its  errors  by  including  the  errors  8rt  and  8ra  in 
Equation  (86)  to  obtain 


r  =  'n/  [Pn  +  Sr"  -  Sr"]T  [pn  +  Sr"  -  8r“] 

Rearranging  Equation  (119)  and  dropping  second  order  error  terms  yields 


(119) 


2(8r"-8r")V 


;  =  V[(Pn)T  Pn] 


1  + 


[(pn)Tpn] 

Using  the  approximation  (1  +  x)m  *  (1  +  mx)  for  Ixl  «  1,  and  rearranging,  Equation  (120)  reduces  to: 
(pn)T  [Sr"  -  8r"] 


(120) 


r  =  r  +  ' 


V[(pn)Tpn] 

r  =  r  +  (u£)T  [8rtn  -  Sr^] 


(121) 


(122) 


The  ESTS  indicated  unit  line-of-sight  vector,  u  ",  can  be  expressed  as  the  sum  of  the  true  value  and  its  error,  i.e. 


—  n  n  o  n 
u  =Up  +  8Up 


(123) 


Derived  next  is  a  linearized  equation  for  8ujJ.  Substitute  Equation  (118)  into  Equation  (89)  to  obtain 

pn  +  8pn 


—  n  . 

u  = 


P  V[(Pn)T  +  (8pn)T][pn  +  8pn] 

Rearranging  Equation  (124)  and  dropping  second-order  error  terms  yields 


(124) 


—  n 
U  = 


(pn  +  Spn) 


1 


(pn  +  8pn) 


p  V[(Pn)V  +  2(pn)T8pn]  \  /  2(p”)T8p" 

[1+  (pn)TPn 


(125) 


Using  the  approximation  (1  +  x)m  *  (1  +  mx)  for  Ixl  «  1,  Equation  (125)  reduces  to: 

Rearranging  Equation  (126)  and  dropping  second  order  error  terms  we  obtain: 

— n  n  1. 


(126) 


uj  =  uj  +  ^[I-uJ(uJ)T]8p» 
uj  =  +  “[  I  -  (uJ)TK8r"  -  5r”) 

— n  n  -  n 
Up  =  Up  +  8Up 


(127) 


where 


Su£  =  “[  I  -  u£  (u£)T](8r"  -  8r^). 


(128) 
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The  INS  indicated  Earth  referenced  velocity  v",  can  be  expressed  as  a  sum  of  the  true  value  and  the  error  value, 
i.e. 


-n  n 
va  =  va 


+  Sv 


(129) 


Expressed  ?  in  terms  of  its  errors  by  including  the  errors  8ujJ  and  the  INS  indicated  velocity  error,  8va,  in 
Equation  (86)  to  obtain: 

f=-[(u;)T  +  (8u;)T][v"  +  8v^]  (130) 

Rearranging  Equation  (130)  and  dropping  all  second  order  error  terms  yields 


f  =f-[(8uJ)Tv"]-[(u^T8v"] 


Substituting  (8u£)T  from  Equation  (128),  we  obtain 


f  -  ?  -  (v")T  7  [I  -  uJ(uJ)T][8rtn  -  8r"]  -  (i#T8v" 
Equation  (132)  can  be  written  more  compactly  by  defining 

L^[I-u>5T 

to  obtain 


r  *  ?  +  (v")TL8r"  -  (v")TL  Sr"  -  (uj)x8v" 


(131) 


(132) 


(133) 


(134) 


8.3  Linearized  INS  Indicated  Azimuth  and  Elevation  Angles  Measurements 


The  INS  indicated  measurements  of  the  azimuth  angle  A  and  the  elevation  angle  E  are  expressed  as  the  sum  of 
the  true  value  and  the  error. 

A  =  A  +  5 A  (135) 

E  =  E  +  5E  (136) 

The  development  that  follows  next  is  only  for  the  azimuth  measurement;  but  it  also  applies  to  the  elevation 
measurement. 

To  get  an  expression  describing  8A,  form  the  difference  between  cos(A)  and  cos(A)  and  call  it  Acos. 

Acos  =  cos(A)  -  cos(  A)  (137) 

Acos  =  cos( A  +  SA)  -  cos(A)  (138) 

Using  the  approximation  found  in  Equation  (105),  cos(0  +  SO)  =  cos(0)  -  50sin(0)  for  "small"  86 ,  Equation  (137) 
reduces  to: 
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Acos  =  -5Asin(A) 

(139) 

Solving  for  8A  yields, 

-Acos 

5A  “  sin(A) 

(140) 

r  / — n.  t  —  n  ,  n.  t  n  .. 

[(V  UDV-^V  "py] 

6A  “  sin(A) 

(141) 

..  tiFK  oOXv 

-sin(A)  +  -sin(A) 

(142) 

At  this  point  it  becomes  necessary  to  derive  an  expression  for  5u^y. 

and  recognizing  the  fact  that  u£y  =  u£y  to  obtain  the  following 

It  may  be  found  by  subtracting  u£y  from  u£y 

8upy  =  upy-upy 

(143) 

8u;y  =  (^-^)  u£y 

(144) 

8uJy  -  (-^  -  C$BcJ)  ujy 

(145) 

&4-KV)*-CS(p)*<i]i4 

(146) 

8uJy  -  ujy 

(147) 

K-'i*? 

(148) 

T 

where  the  superscript  *  denotes  the  skew-symmetric  form  of  a  vector  and  y  =  [yx,yy,yz]  ,  the  vector 
p  =  [Px,fS  and  {f)*  =  (y/)*  +  c£(P)*cjj].  Substituting  Equation  (148)  and  (128)  into  Equation  (142)  yields 


X'r'y^z 

the  expression 
8A  ~ 


(i#7  (u>7  +  j  (ugv)T[  I  -  <  (u^)T3(8r”  -  8# 


-sin(A) 


(149) 


In  [1],  it  is  shown  that  the  gradient  of  A  in  the  navigation  frame,  denoted  VA",  is  expressed  by  the  following 


-T  (C)T[I-uX)t] 


(VA  )  = 


-rsin(A) 


(150) 


Using  the  result,  Equation  (149)  can  be  written  more  compactly  as 

The  first  term  in  Equation  (151)  may  also  be  expressed  in  terms  of  VA"  after  the  following  manipulations: 


(151) 
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(152) 

(153) 

(154) 

(155) 

(156) 

(157) 

(158) 

(159) 

(160) 


where  •  and  x  denote  the  inner  and  cross  products,  respectively.  Equations  (154)  and  (156)  were  obtained  by 
direct  application  of  the  vector  triple  product  expansion.  Substituting  the  result  of  Equation  (160)  into  Equation 
(151)  yields 

8A  =  -r(VA")T(uJ)*yr-  r(VA”)T(uJ)*c£/3  +  (VA")T  5rtn  -  (VA")T  8r"  (161) 


By  a  method  similar  to  that  described  above  for  the  azimuth  measurement,  the  corresponding  equations  for  the 
elevation  measurements  are  given  by: 

5E  -  -r(V£")T  (uj)*yr-  r(V£")T(uJ)*<i3  +  (V£")T  5rtn  -  (VEn)J  8r"  (162) 

where 


Thus,  the  linearized  INS  indicated  measurements  of  the  azimuth  and  elevation  angles  are  given  by: 


(163) 
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(164) 


A  *  A  -  r(VA”)T(uJ)^-  r(VAn)T(u|J)*c£/J  +  (VA")t  5r"  -  (VA")t  Sr* 

E~E- r(V£")T(uJ)*yr-  r(V£”)T(uj})*c£/?  +  (V£")T 8rtn  -  (V£")T  8r” 
respectively. 


(165) 
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9.0  Kalman  Filter  Measurements 


This  section  derives  the  equations  for  the  Kalman  filter  measurements.  Also,  derived  here  are  the  linearized 
measurement  equations  needed  for  designing  a  Kalman  filter. 


The  inputs  to  the  Kalman  filter  are  the  difference  between  the  INS  indicated  measurements  and  the  radar 


measured  values,  i.e., 

Zj.  =  r  -  r 

(166) 

•  • 

=  r-r 

(167) 

zA  =  A  -  A 

(168) 

zE  =1  -E 

(169) 

In  a  Kalman  filter  design,  it  is  necessary  to  linearize  any  non-linear  measurement  equations.  The  linearized 
Kalman  filter  measurement  may  be  obtained  by  differencing  the  linearized  indicated  measurements  and  the  radar 
measurements  to  obtain: 


*T  =  (‘O’W  - 8r”]  -  SrCL  -  r8C  +  +  vrQ  +  vrT  + 

Zf  =  (v”)T  LSr"  -  (v")T  LSr"  -  (uJ)T[8v"]  -  SrD  -  fSC  -  f8F  +  +  v-q  +  v^  + 

zA  =  -r(VAre)T(uJ)*  i/a-  r(VAn)T(uJ)*c£j3  +  (VA")T  8r"  -  (VAn)T  Sr"  -  8ATCMp 

+  va  +  V.  +  VA  +  VA  +  VA  +  V.  +  VA  +  VA 

rtMA  aFL  aPU  rtCC  rtMS  ^TN  rtQ  rtRAD 

zE  =  -r(V£")T(uJ)*r-  r(V£”)T(u;)*0  +  (V£")T  Sr"  -  (V£")x  8r"  -  8ETCMp 

+  Vp  +  Vp  +  Vc  +  Vp  +  Vp  +  Vp  +  Vp  +  Vp 

CMA  “fl  'i'U  cCC  ^MS  “TN  eQ  ^RAD 


(170) 

(171) 


(172) 


(173) 
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10.0  Error  State  Differential  Equations 

The  error  state  differential  equations  for  an  INS  have  been  presented  in  [4,  13,  14,  15].  The  result  of  this  work  as 
well  as  the  results  presented  in  this  report  concerning  radar  error  states  are  presented  in  this  section. 


The  basic  definition  for  the  error  state  vectors  associated  with  the  INS  indicated  outputs  and  the  SAR  and 
monopulse  radar  measurements  are  given  as: 


*lnd  =  [(5r")T((5v")T/VTkvT  eT  k/  (Sr")T  pT  ]T 

(174) 

xRadar  =  t&CL  ^ATEMP  ^TEMP^ 

(175) 

where  V  is  the  accelerometer  bias  vector,  ky  is  the  accelerometer  scale  factor  bias  vector,  e  is  the  gyro  bias  vector, 
and  kg  is  the  gyro  scale  factor  bias  vector. 


Generally,  the  error  states  are  modeled  as  linear  differential  equations  of  the  form: 

x(t)  =  F(t)x(t)  +  w(t)  (176) 

T  T  T 

Consider  the  augmented  state  vector  x  =  [xjnd  .tRadar]  considering  of  both  the  indicated  output  error  states  and 
radar  error  states.  Using  this  augmented  vector,  Equation  (176)  may  be  written  as: 


- 1 

HT1 

CX 

i 

3 

o 

_ i 

xInd 

r  WInd 

L  xRadar  _ 

a? 

o 

_ 1 

_  xRadar  _ 

_  wRadar  . 

(177) 


The  fundamental  dynamic  matrix  FInd’  is  given  in  Table  1. 


Table  1.  Indicated  output  error  state  dynamic  matrix  FInd~ 
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0 
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0 

• 

8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

K 

0 

0 

0 

0 
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0 

0 

0 

o 

0 
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0 

0 
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0 

0 

o 

0 

0 

0 

0 

0 

where, 


F  = 


2 

CO2  -  oos  -  D2  -  (cox  +  nx)px  px  -  [coy  +  Qy]  p 

o>2-cOs-Q2-(co  +Qv)p. 


'Pz  ‘  ^x]  Py 

Pv  *  [®x  "*■  Pz 


’s  \wy.-.y/Kx 
-p  -  [cov  +  QJ  pz 


-py-[(oz  +  Qz]Px 

Px-I®z  +  £2z]Py 
(o2-cOs-£22-(co l+Q  )p 


(178) 


F,„,  = 


L  Vlf/  ‘ 


(0  = 


o 

1 _ 

2®z 

-2toz 

0 

2(oy 

-2©x 

r  o 

Az 

> 

0 

L  Ay 

A 

0 

-coz  1 

CO 

z 

0  - 

£ 

CO 

X 

-2coy 

2co„ 


-®x 

0 


(179) 


(180) 


(181) 


where  A  is  the  acceleration,  p  is  the  angular  rate  from  the  Earth  fixed  coordinates  to  the  navigation  coordinates,  Q 
is  the  Earth's  angular  rate  relative  to  inertial  coordinates,  (0  is  the  angular  rate  from  inertial  to  platform 

coordinates  (to  =  Q  +  p),  (i\  =  g/R  is  the  Schuler  frequency  squared,  g  is  the  gravitational  acceleration,  R  is  the 

Earth's  radius,  oo2  =  coTco.  The  subscripts  x,  y,  and  z  correspond  to  the  directions  in  the  navigation  frame.  Finally, 
the  d  superscript  on  a  vector  is  used  to  denote  a  diagonal  matrix  where  the  diagonal  elements  consists  of  elements 
of  the  vector,  e.g. 


Aa  = 


Ax  0 
0 


0 

0 


0  0  A, 


(182) 


The  matrix  FRadar  for  the  radar  errors  is  given  in  Table  2.  The  vectors  used  to  form  the  noise  vector  w  are  defined 
below: 
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Table  2.  Radar  error  state  dynamic  matrix  FRadar. 


^rCL 

8fD 

SC 

SF 

^  TEMP 

®  TEMP 

^rCL 

0 

0 

0 

0 

0 

0 

8rD 

0 

0 

0 

0 

0 

0 

8C 

0 

0 

0 

0 

0 

0 

8F 

0 

0 

0 

0 

0 

0 

^atemp 

0 

0 

0 

0 

^Atemp 

0 

• 

^TEMP 

0 

0 

0 

0 

0 

^Etemp 

wlnd  =  [°T  WV  WP  0T  0T  0T  0T  0T)T  (1 83) 

*W.=  [»00OVVlT  (184) 

where  wv  is  the  accelerometer  white  noise  vector,  w^  is  the  gyro  white  noise  vector. 
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11.0  Linearized  Measurement  Model  Summary 

In  general,  a  linearized  measurement  model  is  given  by  equations  of  the  form: 

z  =  tfx  +  v.  (185) 

where  H  is  the  measurement  matrix  and  v  is  the  measurement  white  noise.  This  section  summarizes  the 
linearized  measurement  model  derived  in  this  report  by  expressing  the  results  in  the  form  of  Equation  (185). 


Z  =  [  ^Ind  ^Radar  ] 


XInd 

_  xRadar  _ 


+  V 


(186) 


Tables  3  and  4  provide  a  summary  of  the  error  measurement  matrices  #Ind  and  //Radar.  The  measurement  white 
noise  vector  v  is  given  by: 


1  +  v  +  v  +  v 

rDES  *Q  lT  *OT 


+  v;  +  v:  +  v; 


DES 


4OT 


VA  +  VA  +  VA  +  VA  +  VA  +  V 


+  VA  +  V  A 


■MA  nFL  rtPU  nCC  rtMS  AQ 


Vp  +  Vp  +Vp  +  Vp  +  V! 
^MA  nFL  ^PU  nCC 


P  -r  ^p  +  Vp  +  Vp  +  Vp  | 

Erv'  nMS  IZTN  ^RAD  -J 


(187) 


Table  3.  Error  Measurement  matrix  (//Ind). 


12.0  Summary  of  Errors  Budget 

Table  5  contains  a  proposed  error  budget  for  an  INS  and  SAR  system.  This  table  was  compiled  from  data  available  it 
[16]  for  the  Honeywell  H423  strapdown  inertial  navigation  system.  Tables  6  and  7  contain  a  summary  of  the  erroi 
budget,  described  in  this  report,  for  SAR  and  monopulse  radar  measurements,  respectively.  These  tables  were  compilec 
partially  from  a  composition  of  data  from  the  available  literature  and  partially  from  experience  and  engineering 
intuition. 


Table  5.  INS  Error  Budget  Summary- 


Error 

Error  Description 

RMS  Value 

Correlation  Time 

V 

Acceleration  Bias 

30  Hg  . 

bias 

ky 

Accelerometer  Scale  Factor  Bias 

175  PPM 

bias 

Wv 

Accelerometer  White  Noise 

0.0025  FPSVHz 

white 

e 

Gyro  Bias 

0.004  Deg  /Hr 

bias 

ke 

Gyro  Scale  Factor  Bias 

2  PPM 

bias 

"V 

Gyro  White  Noise 

0.0020Deg/Hr/\/Hz 

white 

Table  6.  SAR  Geometry  Error  Budget  Summary. 


1  Error  | 

Error  Description 

RMS  Value 

Correlation  Time 

Range  Clock  Error 

0.01  feet 

bias 

8?d 

Doppler  Measurement  Error 

0.001  feet/sec 

bias 

!  8C 

Propagation  Speed  Error 

10  PPM 

bias 

SF 

Frequency  Error 

20  PPM 

bias 

vr 

rDES 

Range  Designation  Error 

2  Ar 

white 

Vs 

Range  Rate  Designation  Error 

2  At 

white 

% 

Range  Quantization  Error 

Ar/yfn 

white 

Vio 

Range  Rate  Quantization  Error 

Ar/yjn 

white 

\ 

Range  Timing  Error 

V 

white 

\ 

Range  Rate  Timing  Error 

V 

white 

V 

Range  Miscellaneous 

0.5  Ar 

white 

V 

Range  Rate  Miscellaneous  Errors 

0.5  Af 

white 
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Table  7.  Monopulse  Radar  Error  Budget  Summary. 


Error  | 

Error  Description 

RMS  Value 

Correlation  Time 

■Si! 

Antenna  Roll,  Pitch,  and  Yaw  Errors 

1.20  mrad 

bias 

^atemp 

Azimuth  Temperature  Induced  Errors 

0.40  mrad 

bias 

®TEMP 

Elevation  Temperature  Induced  Errors 

0.40  mrad 

bias 

i 

Azimuth  Mechanical  Alignment  Error 

0.30  mrad 

white 

> 

Elevation  Mechanical  Alignment  Error 

0.30  mrad 

white 

Azimuth  Flexure  Error 

0.10  mrad 

white 

Elevation  Flexure  Error 

0.10  mrad 

white 

Azimuth  Radome  Phase  Uniformity  Error 

0.35  mrad 

white 

%u 

Elevation  Radome  Phase  Uniformity  Error 

0.35  mrad 

white 

VAcc 

Azimuth  Rx  Channel  Cross  Coupling  Error 

0.15  mrad 

white 

1 

Elevation  Rx  Channel  Cross  Coupling  Error 

0.15  mrad 

white 

m 

Azimuth  Monopulse  Slope  Error 

0.10  mrad 

white 

VEms 

Elevation  Monopulse  Slope  Error 

0.10  mrad 

white 

Azimuth  Thermal  Noise 

0.20  mrad 

white 

Elevation  Thermal  Noise 

0.20  mrad 

white 

X 

Azimuth  Quantization  Error 

0.20  mrad 

white 

1 

Elevation  Quantization  Error 

0.20  mrad 

white 

Azimuth  Radome  Error 

0.40  mrad 

white 

V]Erad _ 1 

Elevation  Radome  Error 

0.40  mrad 

white 
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13.0  Conclusions 

The  main  objectives  of  this  report  were  to  provide  some  basic  background  on  SAR  targeting,  describe  the  SA1 
measurements  and  measurement  errors,  and  present  a  derivation  of  the  equations  needed  to  integrate  navigation  an 
SAR  measurements  in  a  Kalman  filter.  The  focus  of  this  report  is  on  the  integration  of  SAR  with  an  inertial  navigatioi 
system. 

In  this  report,  we  have  presented  an  introduction  to  ’’unfocused  line-by-line”  processed  SAR,  "focused  line-by-line 
processed  SAR,  and  Doppler  processed  SAR.  We  have  explained  how  each  SAR  measurement  is  obtained  and  ho\ 
various  error  sources  affect  the  measurements.  Both  the  SAR  measurement  errors  and  inertial  navigation  errors  wer 
characterized  by  Stochastic  models.  Finally,  we  have  derived  both  the  non-linear  and  linear  measurement  equation  use< 
in  a  Kalman  filter. 

The  proposed  Kalman  Filter  design  provides  an  optimal  means  to  integrate  SAR  and  INS  data.  It  is  capable  of  utilizing 
the  entire  history  of  SAR  measurements  over  multiple  SAR  images  to  improve  both  SAR  and  INS  measurements. 

The  location  of  the  desired  targets  may  or  may  not  be  known  a  priori.  If  the  target's  location  is  known  a  priori  vi 
intelligence  information,  then  the  Kalman  filter's  estimate  is  initialized  with  this  data.  Subsequent  SAR  measuremen 
are  then  used  to  improve  and  refine  the  a  priori  intelligence  data.  When  the  target's  location  is  not  known  a  priori ,  th 
Kalman  filter  is  initialized  with  the  location  computed  from  the  first  set  of  SAR  measurements.  This  method  is  usefu 
when  the  mission  requires  locating  "targets  of  opportunity". 
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